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We use an off-lattice microscopic model for solutions of equilibrium polymers (EP) in a lamellar 
shear flow generated by means of a self-consistent external field between parallel hard walls. 

The individual conformations of the chains are found to elongate in flow direction and shrink 
perpendicular to it while the average polymer length decreases with increasing shear rate. 

The Molecular Weight Distribution of the chain lengths retains largely its exponential form in 
dense solutions whereas in dilute solutions it changes from a power-exponential Schwartz distribution 
to a purely exponential one upon an increase of the shear rate. 

With growing shear rate the system becomes increasingly inhomogeneous so that a characteristic 
variation of the total monomer density, the diffusion coefficient, and the center-of-mass distribution 
of polymer chains of different contour length with the velocity of flow is observed. 

At higher temperature, as the average chain length decreases significantly, the system is shown 
to undergo an order-disorder transition into a state of nematic liquid crystalline order with an easy 
direction parallel to the hard walls. The influence of shear flow on this state is briefly examined. 
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1. Introduction 



Systems in which polymerization takes place under condition of chemical equilibrium between polymer chains and 
their respective monomers are termed "equilibrium polymers" (EP)tH. The interest to EP from the point of viewj-pi 
both applications and basic research has recently triggered numerous investigation, including computer simulationsO'Ej 
in an effort to avoid difficulties with laboratory experimentsu and approximations as the Mean Field Approximation 
(MFA). Recently the bask scaling concepts of polymer physics were tested by extensive Monte Carlo (MC) simulations 
of flexible EP on a latticea The results suggest that despite polydispersity, EP resemble conventional polymers (where 
the polymerization reaction has been deliberately terminated) in many aspects. However, dynamic aspects of their 
behavior may still be very different: for example, the constant process of scission and recombination in EP gffers an 

under wayu. 



Recently observed 
IB are not completely 



additional mechanism of stress relaxation!]. Computer experiments on EP dynamics are alrea* 

Considerably fewer simulation studies of non- equilibrium properties of EP have been report 
phenomena such as shear banding structure, shear inducing structure and phase transitions! 

understood. An earlieptheoretical workt£l, for instance, predicted a decrease in average size of dilute rod-like micelles 
whereas a later studyEZl concluded that rod-like micelles should ppw at higher shear rates. Since it is known that 
viscoelastic surfactant solutions show unusual nonlinear rheologvES, it is clear that much more research in this field 
is needed before complete understanding of nonlinear properties of EP is achieved. 

Since EP behave in many respects as conventional "dead" polymersu, comparisons with the latter where much 
more work on shear flow effects has been done so far, could prove very useful. Thus inhomogeneity of flows, due 
to the presence of boundaries, and its impact on polymer behavior may be directly observed experimentally by 
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means of evanescent wave-induced fluorescence methodliij that can probe the polymer concentration in the depletion 
layer adjacent to the walls. Coil stretching of dilute flexible ppljaaers in a flow, diffusion and density profiles as 
well as "slip" effects near walls have been treated theoreticallyc3~0 and by computer simulationsc-3~E!3, and as we 
shall demonstrate below, many of these early results compare favorably with what we observe for EP in the present 
investigation. 

In the present study we employ a dynamic Monte Carlo algorithm in order to study EP properties in shear rate. 
The flow of the system in a semi-infinite slit of thickness D is induced by applying an external field F with magnitude 
which changes linearly across the slit and is parallel to the hard walls of the container. Thus the jump rate of the 
monomers becomes thus biased along the x-axis and a flow of the system through the periodic boundary sets in. 

One should emphasize that such an investigation should focus on the linear response in a laminar shear flow. MC 
methods cannot account for hydrodynamic interactions in principle and the transition from laminar to turbulent flow 
can be simulated by means of Molecular Dynamics (MD) only. The linear response breaks down at field intensities 
when the maximum flow velocity is attained, i.e. when all 100% of the random jumps along the cc-axis are forced 
to occur, say, in positive direction. Any further increase of the field F will then fail to accelerate the particles any 
further. Even with these limitations, however, is appears that this kind of MC simulation of EP in a shear flow is 
warranted, given the considerably longer time periods or systems sizes a MC methods may handle as compared to 



All Monte Carlo studies of EP so far have been performed on suppbic lattice either exploiting an analogy of the P<rttj= 
model of magnetism to random self-avoiding walks on a latticetdO, or using the Bond Fluctuation (BFL) Modelcra. 
These lattice models were developed and extensively used for monodisperse systems of conventional polymers and are 
known to faithfully reproduce their dynamic (Rouse) behavior. For the purpose of shear flow studies a disadvantage 
of these models , due to the discrete structure of the lattice, appears obvious: monomers would block each other on 
the lattice at higher shear rates. Random jumps would have to be of the size of single monomers only, and, last not. 
least, the artificial cubic symmetry would predetermine ordering effects along the three major axes of the latticeEZI 
thereby questioning possible phase transitions into liquid crystalline order. 

In the present work we employ an off-lattice model of EP, designed to overcome these and other shortcomings 
of previous lattice models and to serve in examining the role of polymers (semi)-flexibility. An off-lattice model 
should be a better tool in dynamic studies of a broader class of soft condensed matter systems where bifunctionality 
of the chemical bonds might be extended to polyfunctional bonds, as this is the case in gels and membranes. A 
comprehensive comparison of this off-lattice algorithm to earlier lattice modelstl shows that all properties of EP 
derived in former investigations, are faithfully reproduced in the continuum too. 



2. Description of the model 

As in our earlier off-lattice bead-spring model of conventional polymer chainSEl, a coarse-grained polymer chain 
consists of I beads or " effective monomers" . These are connected by springs which represent " effective bonds" and 
are described by a FENE (finitely extendible non-linear elastic) potential: 



Ufene{t) = oo, otherwise 

where r is the distance between two successive beads, ro = 0.7 is the unperturbed bond length with maximal extension 
Imax, R — Imax — To = 0.3, and fc/2 = 20 (in our units of energy fc^T = 1.0) is the elastic constant of the FENE 
potential which behaves as a harmonic potential for r — ro <C R- Thus Ufene{t ~ To) ~ — § ( r ~ r o) 2 but diverges 
logarithmically both for r — > l m ax and r — > l m in = 2r — l m ax- We choose our unit of length such that l max — 1 and 
then the hard core diameter of the beads l m in — 0.4. All lengths as, e.g. the linear size of the simulational box, arc 
then measured in units of l max . 

According to Eq. ([!]) the net gain of energy of a monomer which forms a bond with a nearest neighbor at distance 
ro is then equal to the "bond" energy J. In EP these strong attractive bonds between nearest neighbors along 
the backbone of a chain are constantly subject to scission and recombination. In the present model only bonds, 
stretched a distance r beyond some threshold value, r^reak = 0-8l ma x, attempt to break so that eventually an energy 
Ufene{t) > in the interval between and J could be released if the bond is broken. 



MD. 




for 




(1) 



2 



Each monomer has two unsaturated bonds which may be either engaged in forming a strong saturated bond 
between nearest neighbors along the backbone of a chain (when the originally unsaturated bonds of such neighbors 
meet and become a parallel pair) or remain free (or " dangling" ) as in the case of chain ends or non-bonded single 
monomers. In order to create a bond, however, the respective monomers must approach [j] each other within the 
same interval of distances ri, rea k < r < l max where scissions take place. While covalent bonds are thus constantly 
broken or created during the simulation, we would like to emphasize that no formation of ring polymers is allowed 
and this condition has to be observed whenever an act of polymerization takes place. 

The non-bonded interaction between monomers is described by a Morse potential, 



U M (r) = exp [-2a (r - r min )} - 2 exp [-a (r - r mm )] , for < r - r min < oo, 



(2) 



where r mm = 0.8, and the large value of a = 24 makes interactions vanish at distances larger than unity, so that 
an efficient link-cell algorithms for short-range interactions can be implemented. In the present study we maintain 
our system in the "good solvent" regime and, therefore, only keep the repulsive branch of Eq.(||), shifting it up the 
positive y-axis so that Um{t) = for r > r min . The radii of the beads and the interactions, Eqs. ([!]), (0), have been 
chosen such that the chains may not intersect themselves or each other in the course of their movement within the 
box, so that 'excluded- volume' interactions as well as the topological connectivity of the macromolecules are allowed 
for. 

We introduce the shear rate B by defining an external field F whose only component is directed along the x-axis 
and changes linearly along the z-dimension of the box - Fig^: 



F x (z) = B(z - Z max /2), 



dF x (z) 
dz 



= B 



(3) 



so that the bias changes sign at the middle of the box Z max /2. 



Fig. 1. Definition sketch for external field variation, F x (z), causing flow in x-direction between infinite parallel plates. 



A standard Metropolis algorithm governs monomer displacements, whereby an attempted move of a randomly 
selected particle in a random direction is taken from a uniform distribution within the interval — ^ < Ax, Ay, Az < i. 
The presence of impenetrable walls at z — and z — Z max is observed by rejection of all those jumps of the monomers 
which would otherwise cause them to leave the box through the planes at the bottom and the top. Thus jumps are 
attempted with probability 



1, 
0. 



for 



\ < Ax < | 
otherwise 



Patt(Ax) -- 

and accepted with probability, equal to 



(4) 



1 Recombination for r < rbreak would violate detailed balance if scissions occur at r > r break only. 
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(Ax) = I ex Pi-( E new + AW - E old )/k B T], for E new + AW > E old , _ 
1-1 " ' 1, otherwise, ^ ' 

where E new and E a id are the energies of the new and old system configurations, and AW is the work, performed by 
the field when a monomer jumps from position r « to a new position, r new : 

AW = / Fdr = / F(z(x))dx. (6) 
J fold Jo 

Since the potential of the field F, Eq.@, is not scalar, AW depends on the path followed by the monomer during 
a jump r id — > r new . During a jump this path z(x) is a straight line, = Az x _|_ ^ oW) so that from Eq.(||) one 
obtains for the work 



AW 



[ AX F x .dx = B Znew + Zold Ax = FAx, (7) 
Jo 2 



where F denotes the average from the values of the field in positions f id and r new . With F, as defined in Eq.|^, 
one satisfies the condition of microscopic reversibility with respect to the movements of the particles. Note that it 
is the microscopic reversibility which requires the use of F, rather than F a id, for instance, in the determination of 
AW. This becomes immediately obvious by considering the displacement of a single non-interacting particle in the 
external field: for a self-consistent algorithm, the update rules for the forward and reverse moves have to be identical 
(time reversible). 

The field, F, is thus introduced in the system as an additional term in the Boltzmann probability in Eq.(^|) whereby 
the probability for jumps along or against the field becomes then strongly biased. This term makes the energy in the 
Boltzmann factor in Eq.(||) a monotonously decreasing function of x. The Metropolis algorithm tends to find the 
minimum of this energy thus continuously driving the system to an unreached minimum which gives a continuous 
flow in x direction. One can readily estimate the average jump distance, Sx, if AW is assumed to be larger than the 
microscopic interactions Ufene and Um- Setting E new w E a id = for simplicity in Eq. (|5|), with F « F one has 

Sx = JxP acc (x)dx = j}xexp(-Fx)dx + f\xdx ^ -1 + [l - (l + f ) Cxp (-|:)]/F 2 
JP acc (x)dx fi ex-p(-Fx)dx + f\dx I + I 1 - cx P(~f )\/ F 

which yields Sx = for F = and Sx sa — | for sufficiently strong fields F, (the sign of Sx depending on which half 
of the box is considered). Thus Sx remains bound from above, no matter how strong the applied bias is chosen^]. 

In order to provide for the reversibility of scission and recombination events, a Monte Carlo time step (MCS) has 
been performed after N particles of the system are randomly chosen and attempted to move at random whereby 
existing bonds are kept as they are. After that N monomers are again picked at random, i.e. one of the two bonds 
of each is randomly chosen and, depending on whether this bond exists and points to an existing neighbor along the 
same chain, or it doesn't, the bond is attempted to break or to create. 

Thus, within an elementary time step (MCS) N random jumps and as many attempts of bond scis- 
sion/recombination are carried out, each subject to the Boltzmann probability that the respective attempt is suc- 
cessful. It is clear that in a system of EP where scission and recombination of bonds are constantly taking place 
the particular scheme of bookkeeping is no trivial matter. Since the identity of a particular chain, or monomer 
affiliation, is in principle preserved for no more than one MCS, the,-data structure of the chains can only be based 
on the individual monomers (or, rather, links) as suggested recentlyu. Thus each monomer has two links. The links 
associated with a given monomer are pointers which may either point In the latter case a link then represents an 
unsaturated dangling bond. Thus a large number of particles may be simulated at very modest operational memory. 
Results in the present study involve systems of up to N = 65536. The simulational box (slit) sizes are typically 



2 Of course, since the microscopic interactions are always present, this result remains a rough estimate only whose validity 
can be checked by the simulation. 
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16 x 16 x D where D is the width of the slit in ^-dimension. The total density of the monomers c is then defined as 
the number of monomers per unit volume. 

During the simulation the whole system is periodically examined, the number of chains with chain length I, square 
end-to-end distance, i?g, gyration radius, R 2 g , center-of-mass coordinates, displacements, etc., are counted and stored. 
Because of the semi-periodic (in x- and y-direction) boundary conditions the interactions between monomers follow 
the minimum image convention. The computation of the conformational properties of the chains as R 2 , for instance, 
then implies a restoration of the absolute monomer coordinates from the periodic ones for each repeating unit of the 
chain. 

Technical details of this new algorithm will be presented elsewhere, here we will note only that the high efficiency 
in code performance is achieved by extensive implementation of integer arithmetic in this off-lattice model based 
largely on binary operations with variables. Thus, for example, the most heavily involved (modulo) operations which 
provide periodicity of coordinates and the minimum image computation for distances turn out to be redundant. 



3. Simulational Results 



A. Velocity Profiles 



As mentioned in Section El we create a shear flow in our system by applying an external field with constant gradient 
along the z-axis of the box so that the flow is oriented along the x-direction - cf. Eq.(||), parallel to the hard 
walls at z = and z — D. The lower half of the box would then flow in positive, the upper one - in negative 
x-direction. It is expected that in the immediate vicinity of the walls the flow might be somewhat distorted due to 
walls impenetrability. Below, in Fig.(|^)a we plot the mean jump distance per MCS, Sx, measured along the cc-axis 
for different values of B. The z-coordinates of these successful jumps are taken from the respective z-coordinate of 
the monomers. 
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Fig. 2. (a) Variation of the average distance of accepted jumps in x-direction vs z coordinate with the external field 
amplitude B in a box of size Z max = 32. (b) Variation of the average velocity (distance in a;-direction traveled by a monomer 
after 1024 MCS) at density c = 1.0, Z max = 16, and B = 0.1. 



Evidently, in a wide channel with D = 32 only for sufficiently weak field B < 0.3 the average jump distance grows 
linearly with respect to the half- width of the box (for B = it is zero). For B > 0.3 distortions in the 5x profile 
set in because the maximal jump distance is limited to 0.5, as mentioned in Section |[ For a more narrow slit of 
width D = 16 the region of linear response would then extend to higher values of F < 0.7. Therefore most of the 
simulational results in what follows are derived for D = 16. Fig. ||b then demonstrates that the velocity changes 
linearly across the slit for sufficiently small values of the field B. 

In the broad channel, D = 32, at z = 0, D for B = 0.5 one gets F — 8 from Eq.(|3|) so that the average jump 
distance there according to Eq.(||) should be Sx « ±0.178. The value of 5x at the borders of the box, as seen from 
Fig.(||)a, confirms this estimate demonstrating that the role of the microscopic interactions Um and Ufene is small. 

The presence of the walls is felt in their immediate vicinity and some local distortion of the displacement profile 
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appears increasingly pronounced with growing B although it remains spatially contained in a layer of thickness 
roughly equal to monomers diameter. It is interesting to note that this— small increase of 5x (and, therefore, of 
velocity) immediately at the walls resembles the so called "slip effect"c3 , c3 in simple shear flow of dilute polymer 
solutions in a narrow channel. This slip effect can be explained intuitively by the fact that the polymer molecules 
near the wall align themselves more strongly with the flow than those away from the wall, and are thus transporting 
less flow-wise momentum across the flow than would otherwise be the case. Indeed, in our Monte Carlo model 
the attempted jumps which would otherwise bring the monomers through and beyond the walls of the slit are 
always rejected. Since the molecules cannot penetrate the wall, their concentration is reduced at the wall, so that 
their contributiot^-jto, the viscosity is further diminished. Using a model of dumbbells in parallel wall shear flow, 
one can calculatcEilS-Eil both the nonlinear velocity profile of the suspended solutions as well as the center of mass 
concentration profile between the two walls - we shall see in the next section that the latter is qualitatively reproduced 
by our simulational results. 



B. Effect of Shear Rate on Average Chain Length and Molecular Weight Distribution 

We find that the average chain length of the EP solution, L, decreases steadily with growing B, which is in agreement 
with an earlier MD studyJ - Fig. ||: The mean chain length L at this highest shear rate is about 70% of its value for 
a system at rest. Here we should like to point out the existence of considerable fluctuations in the derived values of 
L for F^O - the statistical error has been reduced at the expense of considerable computational effort. Note that 
the reduced L/Lq mean chain length dependence on shear rate B (Lq is the mean chain length of the solution at 
5 = 0) appears to decrease nearly linearly with B: L/Lq = 1 — 0.35-B which simply follows (cf. Eq. (|To|) below) from 
the exponential dependence of L on interaction J: L/Lq — exp[— (J — bB)/2}/ exp(— J/2) = cxp(&£?/2) w 1 — bB/2 
with the constant b — 0.7 measuring the effective decrease of bond energy J due to shear B. Another interesting 
observation is that the rate of decline is apparently independent of density c, at least for the small values of F 
considered in the present work. 
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Fig. 3. Relative decrease of mean chain length L in lamellar flow field versus B at J/kT — 7, densities c = 1.0 and 
c = 2.0 and D = 16. 

The form for the MWD C(l), that is, the concentrations of chains of contour length I, appears to change qualita- 
tively for equilibrium polymers in dilute solutions. This change is in line with the predictions of the recent scaling 
theory of EPq where we demonstrated that the purely exponential form of the MWD, P{x) oc exp(x), corresponding 
to concentration/chain length regimes in which density correlations are suppressed (typically beyond the semi-dilute 
threshold), is replaced by a 'rounded' Schwartz power-exponential distribution: 

P(x)dx - { eX ^~ x) f X (L»L*) 
( ' ~ 1 Wf) xl exp(- 7 x)dx (L < L*) W 

when correlations are important (typically dilute concentrations). In Eq.(^|) the reduced chain length, x — l/L, is 
taken as ratio of the particular chain length / to the mean chain length L, 7 is the critical exponent of the n — > 
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vector model, (in 3D 7 w 1.165 while its mean field value is jmfa — 1), and L* marks the average chain length 
at the crossover from dilute to semi-dilute concentration, (c — > c*), of EP solutions. The mean chain length L was 
predicted and confirmed to vary with dimensionlcss bond energy J/ksT as 

L oc c a exp(SJ/k B T) (10) 

with exponents ay = &d = 1/(1 + 7) ~ 0.46 in the dilute and a s = 1/2(1 + (7 — V)/(yd — 1) » 0.6, <$d = 1/2 in 
the semi-dilute regime. In Fig. [| we plot C(l) for a system at rest (c = 0.5, B = 0) and at maximum shear rate 
(_B = 0.7) to demonstrate that the form of MWD changes qualitatively 




1CT 6 I ' 1 ' 1 ' „ ,„..„ ,. [ 2 " 4.6 S 



100 200 300 X = l/L 

Fig. 4. (a) Molecular weight distribution in a system of EP at rest, B — 0, and in a flow, B = 0.7, at J/kT — 7 and 
monomer density c = 0.5. (b)The same as in (a) in dimensionless units x = l/L fitted with Eq.(H) with 7 = 1.16. 

when shear is imposed. Thus it is evident from Fig. |] that in the presence of shear the correlations in polymer 
concentration in our dilute system of EP are effectively suppressed and the Molecular Weight Distribution is very 
well reproduced by the simple exponential function expected if a MFA description of the system holds. This finding 
can be understood if one recalls that the imposition of an external field with a shear rate B has a twofold effect on 
the polymers: (i) it effectively reduces the bond strength, J —> J — bB, which makes the polymers shorter, and (ii) 
the shape of the polymer coils is changed towards more rod-like shape with the longest axis oriented along the field. 
Meanwhile it is well known that a system of rods exhibits a Mean-Field-like behaviorB which is here manifested by 
the change in the MWD - Fig.|. 

C. Effect of Shear Rate on Chain Conformations 

As the flow becomes faster, the individual shape of the chain_coils change too. For weak shear rates B, the relative 
distortions (as measured for instance by flow birefwHgenceEi) are essentially proppjjtional to tB where r is the 
largest relaxation time of the unperturbed moleculeEil. Brownian dynamics studiesc3 for Hookean dumbbells in a 
steady shear flow confirm this relative increase of the end-to-end distance with the shear rate both with and without 
hydrodynamic interactions included. A confirmation of these early predictions follows from our simulational results 
too - Fig.|. 
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In Figj^ we plot the difference in gyration radii, typical for a large section of the length distribution and averaged 
over all EP, in two characteristic cases - system at rest and in a flow, 




Fig. 5. Variation of coil size along the field, R gx , and perpendicular to field, R gz , with chain length I and two values of 
B = 0.0, and 0.7. The inset shows the total R g in normal coordinates. Here J/kT — 7 and the density c = 0.5 in a narrow 
slit D = 16. The average coil size measured at zero shear rate is R^ = 50.03 and R? g = 8.325. 

indicating that chain coils in a flow become more extended along the field and compressed parallel to it. The result 
is a total increase of Rg as the whole system starts drifting along field's direction. It is also evident from Fig. |E| 
that this asymmetry in the components of Rg becomes progressively more pronounced as I gets larger. The fact that 
Rg Z is somewhat smaller than R 2 gx even at rest is due to the presence of hard walls at z — 0, Z max which slightly 
deforms the coils in z - direction. We have not given here a plot of the z-dependence of R? g because of the surface 
segregation of chain lengths, induced by the parallel plates. This segregation populates the vicinity of the walls with 
single monomers and very short species. In contrast, the longer chains reside at least a distance R g always from the 
walls (see next Section). Such a distribution of centers of mass with respect to chain length takes place in EP even 
at rest (and is further enhanced by the flow) making the MWD a z-dependent quantity and thus interfering with the 
pure effect of coil stretching under flow. 



D. Density Profiles 



The overall transformations which the system undergoes with increasing shear rate, however, become much more 
explicit if density and diffusion profiles are sampled as function of z. This is shown in Fig.^ where the density is 
normalized to unity (J Q c(z)dz = 1). It is evident from Fig.||a that in the absence of bias when the system is at rest 
the total monomer density is uniformly distributed across the box with a typical depletion immediately at the walls 
(at low concentration of the solution). The walls are avoided by the longer chains because of entropic reasons. When 
the system starts to flow a redistribution of density sets in with increasing bias B whereby for the highest shear 
rate one observes a density maximum centered at the middle where the flo"w__velocity is nearly zero. Qualitatively 
this density profile appears to be similar to analytic and simulational resultsEj for the center of mass concentration 
profile between upper and lower walls, obtained earlier for a single dumbbell in a slitQ 

As the concentration is further increased, one observes the onset of typical oscillations in density profiles in the 
vicinity of the walls, Fig.|^b, c. Such oscillations are typical for polymer solutions confined between, flat plates and 
have been comprehensively studied for conventional polymers by Monte Carlo simulations beforeEffij. The observed 
transitions in the monomer density immediately at the walls from a deficit (depletion layer) at low concentration up 
to an excess (for melts) are governed by a competition between entropic and packing effects. Because of the resultant 



3 In the much simpler modela, however, the density profile is independent of the shear rate B. 
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decrease of configurational entropy it becomes unfavorable for polymers to be near the walls. Chains near the walls, 
on the other hand, suffer collisions with the chains away from the walls, and tend to move closer to the walls. At low 
density the entropic effect dominates, while at high densities packing effects prevail. This is clearly seen in Fig.^,,b,c 
for zero shear. In dilute solutions, as seen from Fig.||a, the increase of shear rate leads to effective broadening of 
the depletion layers, adjacpHt to the walls, which is in agreement with EW1F (evanescentr-wave-induced fluorescence) 
experimental observationsc-l of GM, but at variance with an earlier computer simulationEa. 

This density variation across the slit, caused by the shear rate, appears to depend essentially on the overall 
concentration of the system. At larger shear a density maximum still forms for the slowest layer of flow in the box 
middle at concentration c = 1.0 whereas for very dense systems, Fig.^jc, this density redistribution with shear is 
suppressed. One may thus conclude that effects of shear on the density profiles in a slit depend essentially on the 
free volume in the system which is available for rearrangement of the polymer chains. In the broader slit, Fig.||d, 
where the shear rate gradually diminishes in the vicinity of the wall (adjacent layers flow with nearly equal velocity), 
the density profile gets more complex with two local minima and a sharp increase at the walls. This complex picture 
indicates that monomer density is generally increased in locations of zero flow or steady flow with vanishing shear. 




Density o.015Density 
0.025 - t c ' 




Fig. 6. Distribution of total monomer density between the hard walls in a slit of width Z max = 16 at varying shear rate 
(bias) and J/kT = 7: (a) c = 0.5; (b) c = 1.0; (c) c = 2.0; (d) a broad slit with Z max = 32 and c = 0.5. 

Clearly, the changes in these profiles with B when shear sets on reflect some complex reorganization in the 
polydisperse system of EP whereby the "rapids" of the flow may act differently on chains of different length. 

Additionally, even at rest, the system segregates in the vicinity of the walls for entropic reasonsO into layers 
occupied predominantly by chains of decreasing contour length as one gets closer to the wall. These effects are 
indeed seen in FigJ^, where the average positions of the centers of mass of single monomers and of chains of length 
I = 70 are shown at various strengths of the bias field. The single monomers evidently tend to occupy the immediate 
vicinity of the walls and this tendency is enhanced as the shear increases. The long chains, on the contrary, keep at 
distance « R g from the walls while the system is at rest. For growing B their residence is further narrowed around 
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the " slow" region in the center of the box. In view of Fig.^a one may conclude that the deficit of single monomers 
from this region is more than compensated by accumulation of longer chains. In the wide slit (D = 32), only a 
fraction of the long chains still remains in the middle whereas two new maxima at the walls appear. Evidently, this 
happens in those regions where the shear rate for B > 0.4 (cf. FigJ|) nearly vanishes. 





0.6 V 




Fig. 7. (a) Center of mass distribution of single monomers in a slit with D = 16, J/kT = 7 and c = 0.5 at various values 
of the shear rate (bias), (b) Center of mass distribution of chains with I = 70 under the same conditions as in (a). 

One might expect that other properties of the system, related to density, will also be affected by the shear, as 
for instance, the local diffusion coefficient. In Fig.|| we plot a histogram of mean square displacements (MSQD), 
performed by all those particles which remain in the same z-layer within a MCS. 




Fig. 8. Distribution of MSQD during 1 MCS in a box with Z max = 16 for different shear rates (bias) at J/kT = 7: (a) 
for c = 0.5; (b) for c = 1.0. 

The distribution of MSQD, Fig. ||b, develops from being nearly constant (with two small wings at the depletion 
zones) for zero bias to a well defined broad minimum in the middle of the box as B — > 0.7 whereby as a whole it also 
decreases. Evidently, the diffusion profile across the channel reflects simply the variations of the density distribution. 



E. Nematic Ordering of Short Chains 



Most of the simulational results, discussed in the preceding subsections, have been carried out for a sufficiently strong 
energy bond, J/kT = 7, which is equivalent to a rather low temperature of the system. The average chain length 
at J / kT = 7 thereby varies with density within the interval 40 < L < 70 where the flexibility of the chains ensures 
that their conformations correspond to well shaped polymer coils. It is interesting to check whether a change in the 
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mean size of the chains L in some way affects the reorganization of the polydisperse system under shear flow. If 
one reduces the ratio of bond to thermal energy J/fc^T, as mentioned in the previous section, Eq. jlO|), the average 
contour chain length decreases exponentially fast. In the present study we change J/kT from 7 down to 1 whereby 
L drops from « 40 to 2.5. As shown below this leads to dramatic changes in the EP solution. The profiles along the 
z axis for this case of very short chains are shown in Fig. H. 
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Fig. 9. (a) Distribution of the total monomer density across a slit with D = 16 at various rates of shear (bias), (b) 
Average mean square displacement after 1 MCS. (c) Center-of-mass distribution of chains with I — 1,2,3,4,5,10 and 20 for 
B — 0. Each curve is shifted from the previous one along the y-axis at 0.0005 for better visibility, (d) The same as in (c) but 
for B = 0.7. The width of the slit is Z max = 16, c = 0.5 and J/kT = 1. 

The oscillations in the the total density, Fig. |[ suggest a transition of the system into an ordered state of nematic 
liquid crystal with an easy direction parallel to the walls. Evidently the presence of hard walls acts as an external 
ordering held on the chains. The system is dominated by monomers, dimers and other very short species which 
behave largely like stiff rods aligning themselves parallel to the walls. Indeed, Figs. ||c, [)|d demonstrate that this 
ordering is most pronounced for dimers, trimers and tetramers whereas neither single monomers nor chains with 
length I > 5 participate in the ordering. Thus both single monomers, which lack any anisotropy in shape, and longer 
chains with conformations of coils rather than stiff rods are insensitive to the ordering influence of the walls. 

The influence of growing shear rate on the system with L = 2.5 is similar to that in the case of L ss 40 too. From 
Figs, ^c, ^d one may conclude that the longer chains, which are otherwise uniformly distributed, now pull closer to 
the middle of the slit where the flow velocity is zero. This tendency starts with the 4— and 5-mers already and is very 
clearly seen for the 10-mers (there are only few 20-mers in the system at J/ksT — 1 and their statistics is therefore 
rather poor). This effect more than compensates the developing shallow minimum in single monomer concentration 
in the middle of the box. 
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4. Concluding Remarks 

The present simulational study of the impact of shear flow on EP in a slab reveals a number of interesting features: 

• The average chain length in a system of EP decreases steadily with growing shear rate 

• The polymer coil is gradually stretched along the flow direction as the shear is increased 

• The MWD in a dilute solution of EP changes qualitatively when sufficiently strong shear rate is imposed from 
Schwartz- to a mean-field like exponential distribution function. 

• The shear rate introduces inhomogeneity in the system of EP, confined in a slit: the monomer density, the dif- 
fusion coefficient and the concentration of macromolecules with different lengths develop characteristic profiles 
perpendicular to the walls 

• The width of the depletion layer near the wall for long chains grows with increasing shear rate in agreement 
with recent EWIF studies 

Another interesting phenomenon - an ordering transition in a system dominated by the shorter and stiffer chains is 
found to take place upon heating of the system of EP with the result that chain length is reduced. In this case shear 
flow is observed to enhance the degree of ordering in the system. 

Our observations show that the relaxation of a system of EP from a state of rest to that of steady state flow 
is a slow process which requires long time intervals of investigation, probably rendering Monte Carlo simulational 
methods probably more appropriate than Molecular Dynamics. 

We should like to note, however, that the shear rates studied in the present work are limited to low and moderate 
values since stochastic jumps along and against the external field may be biased by means of the Boltzmann factor 
in a MC procedure within the framework of 100% at most. We therefore expect that at even higher shear rates 
the influence of flow on EP properties might be more dramatic. Clearly, additional work and adequate alternative 
methods are still needed to reach comprehensive understanding of the problem. 
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